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The effects of microvascular branching on blood flow were investigated in vivo by microscopic 
particle imaging velocimetry (micro-PIV). We use micro-PtV to measure blood flow by tracking 
red blood cells (RBC) as the moving particles. Velocity flow fields, including flow pulsatility, were 
analyzed for the first four branching orders of capillaries, postcapillary venules and small veins 
of the microvascular network within the developing avian yolksac at embryonic day 5 (E5). 
Increasing volumetric flowrates were obtained from parabolic laminar flow profiles as a function 
of increasing vessel diameter and branching order. Maximum flow velocities increased 
approximately twenty-fold as the function of increasing vessel diameter and branching order 
compared to flow velocities of 100 - 150 micron/sec in the capillaries. Results from our study 
will be useful for the increased understanding of blood flow within anastomotic, heterogeneous 
microvascular networks. 
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ABSTRACT 

The effects of microvascwlar branching on blood flow were investigated in vivo by microscopic parfide imaging wefocimeiry (mtoro-HV), We use micro-PEV lo measure 
blood flow by tracking red blood cells (RBC) as the moving particles Velocity flow fields, including flow pui&sfiliiy. were analysed ior (ho firsl four branching orders of 
capillaries, postcap U la ry venules and small veins of (he microvascular network wilbin (he developing avian yofksac at embryonic day 5 (E5). Increasing volumetric 
flowrates were obtained from the How prdUes as a function of increasing ve&sfll diameter and branching order. Maximum How velocities increased appraWriatety 
Iwehty-fciW will! increasing vessel diameter and branching order relative to flow velocities of 1 00 - ISO pm/sec in Ihe capillaries. Results Iran our study will be useful for 
advancing the understanding of btood flow wilhfrt anastomotic, heterogeneous microvasculaf networks. (Supported by NASA Glenn Ftesearch Center IRD04-54,) 

INTRODUCTION 

If is hypothesized that the dynamics ol blood flow, including velocity, pressure and shear stress, are Influential factors In microvascular remodeling and 
angiogenesis . Blood flow is primarily laminar and therefore by HagemPoiseuiHe's Faw. Wood flow and pressure are proportional to the fourth power of vessel 
radius (r 1 ). Our research focuses on the role of vessel morphology in microvascular remodeling, angiogenesis and lymphangiogenesis. We quantify vascular 
remodeling as a function ol vessel morphology, diameter, branchpoint and branching generation using Ihe NASA Glenn computer software VE5GEN,' * For Ibis 
study, wo are investigating blood flow as a function of vessel diameter by Intravital particle imaging velocimetry (PlV). ■ PlV measures local velocities within a How 
field by frame-to-frame cross-correlation ol particle displacement. We image the optically accessible mlcrovascutature of tire avian embryonic yolk sac. Red blood 
cells (RBC) serve as the moving particles, so that fluorescent tracking particles are not required. 




Figure 1. Microvascular Branclwnt and Inlrjvilnl Imaging. 

(A) Fertilized [’lljiI eggs were (tortured art a'.u in frwSH Pfe? 1 i i 
dishes until emtiryonc day 5 (ESI. Thv yolksac vatfuldtui y. 
which transports bipod purnptd by Ihe embryonic heart is 
cfearty viaiWe on the yolk- {The chprioaHaritoK: membrane 
(CAM) appears as a smal, highly vascularized ha Boon exiting 
Irom Ihe erauryunic gup. (B) A representative venous 
vascular nee or (he ytfk sac comau: ng lour areas selected for 
imaging ol Brl-Br-t <s displayed . 


METHODS 

tntravitat Imaging and Post-Processing, Fertilized quail eggs [Column* column /aporuea) were cultured 
within 6-weli Petri dishes at 37, 5X as described previously.'" At embryonic day 5 (£5), a field of venous 
blood vessels was selected from lire yolk sac vasculature (Fig. 1 A) tor Hie imaging of tour successive vessel 
branchpoints (Sr1-Br4) lhal included successive venous branching, generations and capillary vessels 
(Fig. IB)-. The Fain dish was placed on a microscope stage pre-warmod to 37 U C and also under an infrared 
heal lamp. A digital Image slack (or movie. Fig. 2A) of 375 images (image size of 1024 x 1024 pixels at 125 
frames per sec and 1/300 sec shutter speed) was acquired in grayscale (256 levels) at total magnrficalim of 
144X (resefution = 0.957 prrv'pixel) wain an Olympus SZX12 stereo microscope and a high-speed, 
high-resolulioti Photron Ultima APX-RS CMOS sensor camera and Photron FaslCam Viewer (PFV) 
software, illuminated by g 10GW mercury burner tight source with UV filter. A sot of tour image slacks, 
Brl-Bol, was acquired from each ol five embryonic specimens: one representative image sei is reported 
here. Images were despeckled by subtracting flom each image a spot image obtained by processing the 
image stack wiih slack averaging, baeground sublraclten, and histogram adjustment. Noise was reduced by 
slight Gaussian filtering. A binary vascular mask displaying only Die vascular fidd was generated by reading 
the aligned images as an imageJ stack, calculating (be standard deviation image. and performing limited 
additional noise removal and thresholding (Fig. 26). Images wore aligned using the inverse of Ihe vascular 
mask wilh NIH Imaged pfug-in TurtoeReg. The vascular mask also served lo isdate ihe vaseiiat velocity field 
for analysis by PlV, 



Figure 2. Imaging and Masking cl a Vascular Field, (A) A gray-scale mu-gu Irorn Hit' Br2 image stock 
iBusrrates vaycuiaf irraphctogy and blood flaw (B) Trio binary vascular inask or the flrz image stack was 
used to £ i ) sign via image stack and (2) isolate the vaaoM stem ter analyse of bteod How velocity by PiV. 


Microscopic Parade Imaging Velocimetry (Micro-PIV). Microscopic PlV (micro-PI V) was performed cm Brl -Br4 mage sets using (he NASA PIVPROC software 
program''' lo measure local velocities at each point on Ihe image grid generated by a subregion size of 64x64 pixels and spacing ol 32x32 pixels (lor second pass: 32x32 
pixels and 16x16 pixels, respectively). Background velocities (i.e., velocrbea in tissue outside of Ihe blood vessels) wilhin an image stack were set to zero using the 
vascular mask (Fig. 2B), PIVPROC employs advanced cross-oorrelatjor! and particle backing operations (dflowing background subtraction, and uses luzzy logic tor 
validating (ha identification of correlation peaks and particle pairings h particfe-lradting operations. 

Analysis at Blood Flaw Velocity and Volumetric Flowrate. To resolve flow pulsalilities resulting from periodic velocity cycles due to cardiac pumping of Ihe 
Wood, fhe velocity field gene rated by PIVPROC for each image wilhin a 375-image slack ol Brl-Br4 was averaged over seven lime steps (i.e., three images 
before, and three images attar, the analyte image). The velocity field of blood flow in an image was converted by the scientific visualization software En&ight to a 
colorized map ol velocity vectors. Average velocities were further calculated lor Ihe rise, peak, tall and valley ol the periodic velocity cycle. The volumetric flowrate 
was calculated at each cross-section of a vessel by numerically integrating the velocity profile. 

RESULTS 

Results for a Brt-Br4 image set from a representative experiment with one specimen are illustrated in Figs. 3-4. Cross-sectional velocity profiles in the 
post-capillary collecting veins were approximately parabolic, indicating laminar Wood flow. According to mean values of fhe crossr sections I cuts (red, green 
and blue, Ftg. 3). the volumetric flowrate of 1 .32 pftmin at maximum (peak) velocity In Ihe largest vessel ol Brl (diameter = 169 pm) decreased to 1 .17 pl/min at 
minimum (valley) velocity. The volumetric flowrate at maximum velocity in Ihe largest vessel of Br4 (diameter - 74 pm) was 6.02x1 0 a pl/min. Overall, the 
maximum flow velocity at the center of the vessels for BrFRr4 decreased approximately 2DX. from a maximum of 2CK30 pm/sec for the largest vessel in &rt 
lo 100-150 pm-'scc in smaller capillaries (lor example, see Br4). Absolute differences in the spatial average ol velocity over time declined considerably from 
Brl to Br4 Pulsatile How cycles in Brl - Br4 were approximately equal (0 4 sec) 


DISCUSSION 

By Inlravital micro-Pl V. velocity flow fields were investigated for (he lirsl tour branching orders ol Ihe microvascular network in the quail yolksac at E5. 
Btood flow was measured by tracking RBCs as Ihe moving particles. Parabolic velocity profiles in the larger veins indicate that Wood flow was laminar. By 
measurements ot individual venous vessels, blood velocity and volumetric flowrate decreased with decreasing vessel sue (Fig. 3), which is qualitatively 
consistent with the Hngen-Pniseuilfe law. In future calculations, we will investigate the mathematical relationship between, vessel diameter, velocity and 
volumetric flowrate in our experimental model, By spatial averaging of Ihe velocity flow field (Fig. 3), absolute differences in maximum and minimum 
averaged velocities for periodic flow cycles also decreased wilh decreasing branching order (Fig. 4). We will conclude our study by segmenting capillaries 
and veins in the velocity flow fields of Brl - Br4 using the GRC software VESGEN (Fig. 5).’ * Segmentation (classification) of the vessels into generational 
branching orders will support further morphological analysis of blood flow dynamics as determined by vessel parameters lhal Include vessel diameter, 
branchpoint density, and tortuosity 



Figure 3. Analysis ul Velocity Flaw Fie I its lor Successive Vessel Branching Generations. VotocMy Mow fields irieasuiod by micro- PlV and averaged (or valuers of peak how are illustrated 
tor Srt -Br4 (A-D), where the magnitude of the velocity vectors (|im.’sec) are given by color lookup tobies (LUT) Velocity profitas across the ona|or vassals show good agreement. 
as generated by cr^-sectfootil cute (red. yraan, and blue cuctes). end art? plotted in the upper right corner of the figure panel® The white arrows Indicate me direction or 
venous flow, end the white lines are calculated panicle path fines, tor Hie RBCs. Average velocities; decreased considerably with decreasing vessel diamete r. Velocity profiles are 
highly parabolic, indicating that blood new within the vessels was lamriflr Cyclic puteatilities d Ihe venous blood How drculatJQfl are displayed in graphs or average velocity 
versus lime (gold, center right] 


spatial average of velocity magnitude by frame 
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Figure -i Comparison ol Pulsatile FIdw Magnitude and Cycle far BrlBr4. I ho spatiaJ average ol thu velocity now 
(laid within each masted Image ol Ihe 375 Image slack is plotted versus time In red (Iname-by-liame] and blue 
] 7 - 1 tl s r fmj average] . The averaged velocity magnitude Increased up lo Hvetotii in the flow field containing Ihe 
largest vessels (Br l ], relative- lo the How held ol smallest vessels (Br4). AbsoMe differences in file average 
velocity ol peak-to-vailey blood How decreased ronstoerabty from Bf 1 to Br4. The periodic cycles ol pusatite 
How lor Br1-Br4. however, were aFproKinsalety equal (0,4 sec). 


Figure 5. VESGEN Saltware for Aulomatii; Analysis ot Vessel 
Brandling Generations (G1-G&). I tie VESGEN sottweie analynes 
numerous vessel pisramefere such as vessel diameter, length, 
tortuosity and number density tor each generation ot branching 
within a vascular Iree [soon to be released by NASA Slenn as a 
dowmaadabto sonwaro plug-m lor NJM ImogeJ) ,J As a ima: 
study tor our projod. vetodly How fwids Bri’Br4 measinod by 
micro- PlV will be analyzed according to vessel branching 
generaliort as determined by VE&3EN. 
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